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Abstract 

We investigate the relic abundance of asymmetric Dark Matter particles that were 
in thermal equilibrium in the early universe. The standard analytic calculation of the 
symmetric Dark Matter is generalized to the asymmetric case. We calculate the asym- 
metry required to explain the observed Dark Matter relic abundance as a function of 
the annihilation cross section. We show that introducing an asymmetry always reduces 
the indirect detection signal from WIMP annihilation, although it has a larger annihi- 
lation cross section than symmetric Dark Matter. This opens new possibilities for the 
construction of realistic models of MeV Dark Matter. 
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1 Introduction 



One of the most striking features of cosmology today is that the Universe contains a large 
amount of Dark Matter, whose mass density exceeds that of the known baryonic matter 
by about five times. Cosmic microwave background (CMB) anisotropy observations by the 
Wilkinson Microwave Anisotropy Probe (WMAP) yield an accurate determination of the total 
amount of the Dark Matter p], 

tt DM h 2 = 0.1109 ±0.0056, (1) 

where Odm is the Dark Matter (DM) density in units of the critical density, and h = 0.710 ± 

0. 025 is the Hubble constant in units of 100 km sec -1 Mpc -1 . 

Although we already have a relatively precise measurement of the Dark Matter density, we 
still do not know what the Dark Matter is. Neutral, long-lived or stable weakly interacting 
massive particles (WIMPs) are currently considered as excellent candidates for Dark Matter. 
One of the most popular WIMP Dark Matter candidate is the neutralino in supersymmetric 
theory, which is stable due to the conserved .R-parity [2]. Neutralinos are Majorana particles, 

1. e. they are their own anti-particles. However, this is only one possibility and we do not 
actually have any evidence that the Dark Matter consists of Majorana particles. Indeed, 
most of the known elementary particles are not Majorana particles, but rather have distinct 
anti-particles; moreover, the ordinary matter in the Universe is almost entirely made from 
baryons, with anti-baryons contributing only a tiny fraction of the baryons. It is thus natural 
to consider asymmetric Dark Matter (ADM), for which particles and anti-particles are not 
identical. This allows scenarios where the universe exhibits a Dark Matter asymmetry, i.e. 
there are more Dark Matter particles than anti-particles (or vice versa). 

One motivation for considering this scenario is that the average density of baryons is 
comparable to that of Dark Matter; with Qf, w 0.046, they differ only by a factor of about 
five. This led to speculations that a common mechanism might give rise to the known baryon 
asymmetry and a postulated asymmetry in the Dark Matter sector [3], 0]. 

One attractive feature of the WIMP Dark Matter scenario is that it provides a natural 
explanation for the observed Dark Matter abundance (the 'WIMP miracle"). At early time 
when the interaction rates are high, the WIMP particles are in thermal equilibrium with the 
rest of the cosmic fluid. As the Universe expands and the temperature of the fluid drops, 
the interaction rates become smaller. Eventually the Dark Matter particles decouple from the 
rest of the cosmic fluid, so that their co-moving density is constant ("frozen"). For symmetric 
Dark Matter, the resulting relic abundance is of the same order of magnitude as the observed 
value, if the WIMP annihilation cross section is of weak size. This hints at a connection 
between weak-scale physics (which is also being probed at colliders) and Dark Matter. 
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The abundance of Dark Matter can be calculated by solving the Boltzmann equation 
which describes the time evolution of particle densities in the expanding Universe. In case 
of symmetric Dark Matter, such calculations have been done for various standard [5J |6] and 
nonstandard |S] cosmological scenarios. 

In this work we generalize this calculation to the case of asymmetric Dark Matter. We 
solve the relevant Boltzmann equations both numerically and in an analytic approximation. 
We assume that WIMPs can annihilate only with their anti-particles, so that the number of 
particles minus the number of anti-particles in a co-moving volume is conserved. Without loss 
of generality we assume that there are more particles than anti-particlesj*] We assume that 
the Dark Matter asymmetry, if any, is created well before Dark Matter annihilation reactions 
freeze out. Moreover, we assume standard cosmology during and after WIMP decoupling. 
This implies a constant co-moving entropy density, and absence of late non-thermal WIMP 
production (e.g. from the decay of heavier particles). 

This paper is arranged as follows. In Sec. 2, we discuss the Boltzmann equations and the 
relic abundance of asymmetric Dark Matter. In Sec. 3, we derive an approximate analytical 
formula for the asymmetric case. In general the resulting Dark Matter density depends both 
on the Dark Matter asymmetry and on the Dark Matter annihilation cross section. In Sec. 4, 
we obtain constraints on these parameters from the observed Dark Matter relic abundance. We 
also compute the ratio of anti-particle and particle densities in the allowed part of parameter 
space, and show that a non-vanishing asymmetry leads to a reduced rate of WIMP annihilation 
in galactic haloes. Finally, in Sec. 5 we summarize our results and draw some conclusions. 

2 Relic Abundance of Asymmetric Dark Matter 

Consider a Dark Matter particle denoted by x that is not self-conjugate, i.e. the anti-particle 
X 7^ X- The relic densities of x an d X are determined by solving the Boltzmann equations 
which describe the time evolution of the number densities n x , in the expanding universe. 
Under the assumptions that only XX pahs can annihilate into Standard Model (SM) particles, 
while XX an d XX pahs cannot, the relevant Boltzmann equations are: 



*In other words, we define the "particle" to be the one with the larger density, if an asymmetry exists. 




(2) 




(3) 



2 



where Mpi = 2.4 x 10 18 GeV is the reduced Planck mass, and g* is the effective number of 
the relativistic degrees of freedom. Here we will consider only the case of cold Dark Matter, 
i.e. we assume the Dark Matter particles were already non-relativistic at decoupling. The 
equilibrium number densities n X)Bq and n Xfiq are then given by 



3/2 



where m x is the mass of the WIMP, /i x is the chemical potential of the particles, and g x 
is the number of the internal degrees of freedom of the x particle. We have used the fact 
that fi x = — fi x in equilibrium. As a result, the chemical potential drops out in the product 
n x,eq n x,eq- 111 Eqs-© the term proportional to this product describes XX production from SM 
particles; it is clear that this term should not be affected by the WIMP chemical potentially 

We follow the standard picture of the Dark Matter particle evolution. At high temperature 
the x an d X particles are in thermal equilibrium in the early universe. When T drops below 
the mass m x , the number densities n XjGq and n Xfiq decrease exponentially, as long as m x > |/i x |. 
Eventually the interaction rates V = n x {av) and F = n x (av) therefore drop below H, which 
scales like T 2 , see Eq.Q. The x an d X distributions are then no longer kept in chemical 
equilibrium, and their co-moving number densities approach constants. The temperature at 
which the WIMPs drop out of chemical equilibrium is called the freeze-out temperature. 

The Boltzmann equations (T5]) can be rewritten in terms of the dimensionless quantities 
Y x = n x /s, Y x = n x /s, and x = m x /T, where 

s = (2vr 2 /45)^T 3 (5) 

is the entropy density. If we assume that the universe expands adiabatically during this period, 
the Boltzmann equations become 

dY x X(av) 



dx x 2 

dY x \(av) 
dx x 
where we have introduced the constant 



O^X ^X ^X,eq ^X,eq) ] (6) 
2 O^X ^X — ^X,eq ^X,cq) j v) 



A = —= m x Mpi Jg* 
x/90 X V 



^Here we ignore effects due to Bose enhancement or Fermi suppression in the final state, as usual in the 
treatment of the decoupling of WIMPs. These effects become important only if the Dark Matter particles were 
(semi-)relativistic at decoupling [B]. 
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Subtracting Eq.flH]) from Eq.([7]), we obtain 



This implies 



KXlJU VAiJb 



Y X -Y X = C, (10) 



where C is a constant, i.e. the difference of the co-moving densities of the particles and 
anti-particles is conserved. This follows from our assumption that \ and % only annihilate 
with each other, which could e.g. be due to conservation of some (global) charged Inserting 
Eq. ffHJ]) into Eqs.flS]) and (j7J, our Boltzmann equations become 

f-^« + Cli-?), (12) 

where 

P = Y xm Y xm = (O.U5g x /g*) 2 x 3 e- 2x . (13) 

As noted above, P does not depend on the chemical potential fi x , which simplifies our calcula- 
tion. These equations can be solved numerically; a semi-analytical solution will be presented 
in the next Section. 

In most cases, the WIMP annihilation cross section can be expanded in the relative velocity 
v between the annihilating WIMPg^ 

av = a + bv 2 + C(t» 4 ) . (14) 



If XX annihilation from an S— wave initial state is unsuppressed, the first term in Eq.(p2 
dominates, while for annihilation from a P— wave initial state, a = but in general b ^ 0. In 
all examples we know, annihilation from either the S— or the P— wave (or both) is allowed; 
the expansion (|14p is then sufficient for an accuracy of a few percent. 

In Fig. [1] we show the evolution of the relic abundances for asymmetric Dark Matter \ 
and its anti-particle % as function of x — m/T^ Here we take 6 = for simplicity, with 
a = 5 x 10 _9 GeV~ 2 and C = 10~ n . Results with non-vanishing b are qualitatively the same. 
For comparison, we also show result for the symmetric case (C = 0). We show the actual as 



*C is proportional to the initial (high— T) value of the x ~ X asymmetry, A x = (n x — n x )/(n x + n x ). 
However, A x will change with time or temperature once the WIMPs become non-relativistic, approaching a 
(larger) constant again once the WIMPs decouple. We therefore prefer to use C rather than A x to parameterize 

the x ~ X asymmetry. 

§See ref.[5] for a discussion of scenarios where this expansion does not work. 
^Similar numerical results have been presented in ref, [3], 
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Figure 1: The evolution of the scaled \ an d X abundances as function of x — m/T for a = 5 x 10~ 9 GeV 
6 = 0, m = 100 GeV and C = 10~ n or zero. The \ equilibrium distributions are shown for comparison. 



well as equilibrium values of the x density; for C = (C > 0), the equilibrium \ density is 
equal to (larger than) the equilibrium \ density shown in the figure. Note that x an d X are 
distinct particles even if C = 0. 

As we can see in the figure, at high temperature (small x) Y x and Y% are very close to 
their equilibrium values. This is true irrespective of the initial value Y x (xq) as long as xq is 
somewhat smaller than the freeze-out value xp, xq — xf ^ 1 [8J. Moreover, Y Xfiq (x) is almost 
independent of C as long as Y x>eq (x; C = 0) > C. This inequality is satisfied for 

x<x c ^-HC') + ^\n[-\n(C')], (15) 

with C = 2 5 / 2 7r 7 / 2 ^C/(455f x ) ~ 310C; in the last step we have taken g x = 2, g* = 90. 
Numerically, xc — 24 for C = 10 -11 . xc is therefore also close to the value where the x 
density begins to differ significantly from its equilibrium value, which signals the on-set of 
decoupling. For C = the actual x density then rather quickly approaches a constant. 

For C = 10~ n the x abundance begins to deviate from its equilibrium value also at x ~ 23. 
However, it keeps decreasing for a large range of x— values, only slowly approaching a smaller 
constant value (~ 10~ 13 ) at very large x. The evolution of the x abundance is even more 
remarkable: we have not plotted its equilibrium density in this figure, because for C = 1CT 11 
it almost coincides with the actual x density at all temperatures! The ratio Y x /Y x ^ eq never 
exceeds 1.85; this maximum value is reached for x ~ 25. This can be understood from the 
observation that the x density becomes nearly equal to its equilibrium value, simply given by 
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C, also at very large x, where the \ density is very small. 

This large difference in the evolution of the % and % densities originates again from our 
assumption that only XX pairs can annihilate. Since Y x is bounded from below by C, \ 
particles still find a significant density of partners for annihilation even after the nominal 
decoupling temperature. This explains why Y% keeps decreasing even at rather large x. This 
in turn means that x particles find even fewer partners for annihilation than in the symmetric 
case {C = 0). The x density therefore "decouples" (i.e., T x < H) even earlier than for C = 0, 
although its number density never differs very much from its equilibrium value, as we just 
saw. 




a [GeV 2 ] b [GeV 2 ] 



(a) (b) 

Figure 2: The relic density ilh 2 for particle x an d anti-particle x as a function of the cross section. Here 
we take m x = 100 GeV, g x — 2 and <?» = 90. Panel (a) is for b = 0, while panel (b) is for pure P— wave 
annihilation (a = 0). 



With this understanding of the asymmetric Dark Matter thermal decoupling, we now look 
at how its abundance depends on the annihilation cross section. In Figure [2fa,b) we plot 
the relic density Qh 2 for particle x an d anti-particle % as a function of the annihilation cross 
section, with two different values of C. In panel (a) we assume a constant (S— wave) cross 
section (6 = 0), while panel (b) is for a pure P— wave case (a = 0). 

We see that for small XX annihilation cross section the relic density is independent of the 
asymmetry. This is the case whenever the relic density in the symmetric case (C = 0) is much 
larger than the assumed value of C; in this situation the small asymmetry clearly does not play 
much of a role in the decoupling dynamics. Since for C = the relic density scales essentially 
inversely with the annihilation cross section, this scenario corresponds to small cross sections. 
More precisely, the product of annihilation cross section and asymmetry determines whether 
the asymmetry has significant impact on the final relic density or not. 
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In the regime of small cross section the velocity dependence of the cross section is relevant. 
In the symmetric case the contribution of the b— term in Eq. (I14|) to the final relic density is 
suppressed by a factor xp/3 ~ 7 relative to that of the a— term. A similar scaling holds in 
Fig. [2] for small cross section. 

On the other hand, as already shown in Fig. dj for sufficiently large annihilation cross 
section a non-vanishing asymmetry C does greatly affect the relic density. The scaled x 
density Y x will then simply approach C, while the final x density becomes insignificant. This 
is the limit of extremely asymmetric dark matter. The total DM relic density then becomes 
independent of the value and the velocity dependence of the cross section (although the velocity 
dependence does play a role in deciding what cross sections are "large" in this sense). 



3 Analytical Solution 

While the Boltzmann equations (TTTj) . (Tl2|) can be solved numerically, it is still useful to have 
some analytical solution of these equations, as in the symmetric Dark Matter case. In our 
semi-analytical treatment we closely follow the standard treatment [5] of symmetric Dark 
Matter to solve Eq. fll2p for the x density. The \ density can then be computed trivially using 
Eq.flTO]). 

We begin by introducing the quantity A^ = Y% — Y Xieq . Its evolution follows directly from 
the Boltzmann equation (TT21) : 

^_^_AM [Aj(As + 2}W + GAx] . (16) 

We consider the behavior of the solution of this equation in two regimes. At sufficiently 
high temperature Y% tracks its equilibrium value Y Xfiq very closely. In that regime A^ is small, 
and dA x /dx and A| are negligible. The Boltzmann equation (|T6|) then becomes 

% 5 = "^ (2A Aeq + CA,). (17) 

In order to calculate the left-hand side of Eq.(iT7jl. we need an explicit expression for the 
equilibrium density Y Xfiq (x). By definition the right-hand sides of the Boltzmann equations 
and (J7j) vanish in equilibrium. Hence the right-hand side of Eq. (Tl2|) should vanish as well 
for Y x = Y xm , which implies 

Yl eq + CY^ eq -P = 0. (18) 

This quadratic equation has two solutions, but only one of them yields a positive x equilibrium 
density: 

Y x^ = -^ + J T + P- (19) 



Plugging Eq. fll9p into Eq. (ll7l) and ignoring x relative to x 2 , we have 

2x 2 P 

* ~ X(crv) (C 2 + AP) ■ ( 0) 

This solution will be used to determine the freeze-out temperature for x, which we denote by 

Xp. 

At sufficiently low temperature, i.e. for x > xf, we can ignore the production term oc Y Xfiq 
in the Boltzmann equation (1T61) . so that 

dAy X(av) / a9 ^ a x 

^ = -^( A l + GA x) ■ (21) 

At this point we may assume A x (xp) 3> A x (oo). Integrating Eq. pT]) from to oo then 
yields 

Y x (x ->• oo) = ^— r . (22) 

exp (CX f^(av) x~ 2 dx) — 1 

If we use the non-relativist ic expansion (|14p of the cross section times relative velocity, the 
thermal averaging gives 

(av) = a + 6bx~ l + 0(x~ 2 ) . (23) 



The solution (1221) then becomes 

C 

yi^ oo) = (24) 

x exp [C(47r/V90) m x Mpi ^ (a Xp 1 + 3b x F 2 )] - 1 ' 



where we have used Eq.flE]). Note that for small asymmetry C, or more precisely for small 
argument of the exponential in the denominator, C will drop out, and the x re li c density 
will be inversely proportional to the annihilation cross section. This reproduces the standard 
result |5j. On the other hand, Eqs.(J22J) or (f24|) indicate that the x reuc density will become 
exponentially small if the product of annihilation cross section and asymmetry becomes large. 
We can repeat this exercise to calculate the relic density of x particles. The result is 

C 

Y (x — )• oo) = (25) 

l-exp[-C(An/V90)m x M Pl ^g:(ax F 1 + 3bx F 2 Y ' 



where xp is the inverse scaled freeze-out temperature of %. Note that Eqs. (l24l) and (T25l) are 
only consistent with the constraint ( ITOl) if xp = Xp. Recall from the discussion in the previous 
Section that for sizable C the reaction rate T x is (much) smaller than the rate T x , since x 
particles find (far) fewer partners for annihilation. The fact that our treatment nevertheless 
requires equal "freeze-out" temperatures for x an d X particles indicates that the intuitive 
condition T(xp) ~ H(xp), i.e. reaction rate ~ expansion rate, can no longer be applied 
consistently for the decoupling of asymmetric Dark Matter. 
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For convenience, we express the final abundance in terms of 



tt x h 2 



m x soY x (x —> oo)h 2 



(26) 



Pcrit 



where sq = 2.9 x 10 3 cm -3 is the present entropy density, and p cr it = 3Mpji?o is the present 
critical density. The corresponding prediction for the present relic density for Dark Matter is 
given by 



For small product of annihilation cross section and asymmetry the WIMP mass m x drops 
out in the final result fl27|) . In contrast, if this product becomes large, the x re li c density is 
negligible while the scaled x rei i c number density is simply given by C; in this limit the final 
Dark Matter mass density is thus proportional to the mass of the Dark Matter particle. In 
that regard the Dark Matter mass density of very asymmetric WIMPs behaves like that of 
hot Dark Matter, although our WIMPs remain cold Dark Matter. 

The remaining task is to fix the freeze-out temperature. Here the standard method [5] 
needs some modification. We start from the standard definition of xp which assumes that 
freeze-out occurs when the deviation is of the same order as the equilibrium value of Y x : 



where £ is a numerical constant of order unity. We adopt the usual value [5] £ = \/2 — 1. 
By comparing with numerical solutions, we find that the standard treatment under-predicts 
the x renc density for sizable asymmetry. To very good approximation the ratio of the exact 
analytical solution and the approximation using xp Q as decoupling temperature only depends 
on the product XC(av )(xp )/x Fo , rather than on C, m x and the annihilation cross section 
separately. This is not so surprising, since this expression appears as integrand in Eq.( l2"2"j) . 
Moreover, the ratio of exact and approximate solution remains of order unity as long as the 
above product is < 1. 

We therefore need only a relatively minor modification of the standard treatment. Af- 
ter some trials, we found that the most elegant and simple way to improve the analytical 
approximation is through a small shift of the freeze-out temperature: 



In Fig. [3] we plot the ratio of the exact value of the % particle abundance to our analytical 
approximation. As shown in the figure, the approximate analytic result matches the exact 
numerical result very well as long as XC(av)(xp)/x F < 1. For even larger product of cross 
section and asymmetry the deviation becomes sizable again. However, in that case the x 



n DM h 2 = 2.76 x 10 8 m x [Y x (x ->• oo) + Y x (x ->• oo)] GeV 



(27) 



(28) 




(29) 



9 



m=100GeV,a=10_ GeV_ 

m=30GeV,a=10_ 9 GeV^ 
m=100GeV,a=10 8 GeV 



0.95 



0.95 



m= 1 OGeV , b= 1 GeV_ 
m=3 0GeV,b=10_ 8 GeV_ 
m=100GeV,b=10 8 GeV 



. 9 



0.9 



0.3 



1 . 2 



1.5 



0.1 



0.2 



0.3 



>.*C*b/x F 



Figure 3: The ratio of the exact value of the \ particle abundance to the analytic value of x particle 
abundance. 



relic density is entirely negligible. Quantitatively, our approximation reproduces the exact 
numerical solution to better than 2% (5%) as long as the final x renc density exceeds 10~ 4 
(10~ 9 ) times the final x re hc density. 

Since for large product AC (av) the scaled x re li c number density is simply given by C, our 
simple approximation reproduces the exact x re li c density, and also the exact total WIMP 
(x + X) mass density, to better than 2% for all combinations of parameters. 

We finally note that the shift ( 1291) of xp amounts to < 2% for the parameter range covered 
in Fig. 131 This nevertheless can lead to a significant change of the x re li c density due to the 
exponential dependence of Eq.( l24|) oni f . 

4 Constraints on Parameter Space and Indirect Detec- 
tion Signals 

The Dark Matter abundance can be derived from cosmological observations. This measure- 
ment does not depend on the nature of dark matter, e.g. whether it is symmetric or asym- 
metric, as long as the Dark Matter particles are sufficiently slow ("cold" Dark Matter). Using 
only the CMB data, the WMAP team derived the Dark Matter density given in Eq.(JT]) for the 
minimal ACDM model. The precision could be further improved by including observations of 
baryon acoustic oscillation (BAO) and/or direct measurements of the Hubble constant (Hq) 
into the fits [TJ. We note that besides the statistical error quoted in Eq.([T}), there could also 
be systematic errors, and the result depends somewhat on the model and priors adopted. 
Nevertheless, we may use the measured Dark Matter abundance to place constraints on the 
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asymmetric Dark Matter model. In view of the above remarks, we adopt the conservative 
range 

0.10 < tt DM h 2 < 0.12 (30) 

when setting bounds on the model parameters; this is roughly the ±2cr limit for the abundance 
indicated by the WMAP fit given in Eq.([T]). Note that for asymmetric Dark Matter, the x 
and x contributions have to be added: 



(31) 



In this Section we assume that the annihilation cross section is given by Eq. ffl4l) . We use ex- 
act numerical solutions of the Boltzmann equations; however, the analytical solution described 
in the previous Section reproduces the exact result to better than 5% for all combinations of 
parameters we analyze in this Section. 
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Figure 4: The allowed region in the (a,C) plane for b = (left), and in the (b, C) plane for a — (right), 
when the Dark Matter density fi/i 2 lies between 0.10 and 0.12. Here we take m x = 100 GeV, g x — 2 and 
<7* = 90; the allowed values of C scale to good approximation like l/m x . 



The relation between the cross section parameters a and b (for 6 = and a = 0, respec- 
tively) and the asymmetry C are shown in Fig. [4] for two values of the total Dark Matter 
density. For C = 0, i.e. in the symmetric case, one needs a ~ 4 — 5 x 10~ 9 GeV~ 2 for b = 
(i.e. S— wave dominance), and fe~3x 10~ 8 GeV" 2 for a = (P —wave dominance). These 
cross sections are about two times larger than the cross sections required for self-conjugate 
(Majorana) Dark Matter with g x = 2 degrees of freedom. This is due to the fact that the 
(equal) x an d X contributions have to be added, as in Eq. (l3~il ; each of them is about as large 
as the contribution from Majorana dark matter with the same cross section. The resulting 
doubling of the total Dark Matter density has to be corrected by also (approximately) doubling 
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the annihilation cross sectionjj 

Now consider the asymmetric case, i.e. C ^ 0. For small values of C, the abundance is 
not much affected, so the iso-abundance contours initially rise almost vertically. This can be 
understood from the analytical expressions (|24p and (|25p : if we ignore the shift of xp given by 
(EHD, which is very small for small C, we find Y x [x oo) ~ Y x (x ->■ oo; C = 0) + C/2 + 0(C 2 ) 
and Y x (x — >■ oo) ~ F x (x — >■ oo; C = 0) — C/2 + (9(C 2 ), i.e. the correction to the total Dark 
Matter relic density begins at 0{C 2 ). 

For somewhat larger cross sections the curves quickly flatten out. A larger cross section 
leads to smaller decoupling temperature, and hence to a smaller relic density; this has to 
be compensated by increasing C. Specifically, once the cross section is twice as large as 
that required for vanishing asymmetry, the required asymmetry is already only about 5% 
smaller than the asymptotic value required for large annihilation cross section (i.e., for strongly 
asymmetric Dark Matter, where the \ relic density is negligible and the \ density is simply 
given by C). Since C = gives the minimal allowed value of (av) while (av) — > oo defines 
the maximal allowed value of C, the above statement can be written as 

C((av) > 2(av)\ min ) > 0.95C max . (32) 

Of greater physical interest than the ratio C/C max are the fraction Y x /Y x as well as the 
product aY x Y x , where the particle densities are to be taken at the present epoch, i.e. for very 
large x. These quantities are shown in Figs. [5]and[6l respectively. 

The ratio Y x /Y x determines the contribution of anti-particles to the current Dark Matter 
density. In the left frame of Fig. [5] this ratio is shown as function of m x ■ C, so that the results 
are nearly independent of m x ; the annihilation cross section has been chosen such that the 
relic density has the indicated value (f^M^ 2 = 0.10 [0.12] for the black [grey or red] curves). 
Plotted in this way, the result is almost independent of whether the a— or b— term in the 
annihilation cross section (fl4|) dominates. 

We see that for small C the logarithm of the ratio falls essentially linearly with increasing C, 
with relatively gentle slope. Recall from the discussion of Fig. H]that for constant annihilation 
cross section the ratio would behave approximately like [Y X (C — 0) — C/2]/[Y x (C = 0) + C/2] 
for sufficiently small C; this gives \n{Y x /Y x ) ~ —C/Y X {C = 0). Fig. @] also showed that 
increasing C requires at first only very small increase of the annihilation cross section. 

*One can also compare non self-conjugate Dark Matter with Majorana Dark Matter with g x = 4. Note 
that for fixed annihilation cross section, the final relic density depends on g x only logarithmically (via Xf). 
An asymmetric Dark Matter scenario with given cross section then corresponds to a model of Majorana Dark 
Matter with half the cross section, since in the non self-conjugate case effectively only half the initial states 
contribute to the annihilation (xx an d XX contribute but XX an d XX do not). 
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Figure 5: Ratio of the current anti-particle abundance Y x to the particle abundance Y x as a function of 
m x ■ C (left) or the sum a + 6/6 characterizing the annihilation cross section (right) for different combinations 
of cross sections and total Dark Matter relic density. In the left (right) frame the annihilation cross section 
(asymmetry C) is chosen such that the total dark matter density f^M^ 2 has the indicated value. Results are 
for = 90 and m x = 100 GeV, but are almost independent of m x . 



This behavior continues until C reaches about 95% of its maximal value, where the ratio 
Y x /Y x ~ 0.03. When C is increased even more the ratio plummets. This is caused by 
the required increase of the annihilation cross section, which has to become very large as 
C — > C max . We saw in Fig. |2]that this strongly suppresses the x reuc density. The logarithmic 
slope of the ratio therefore approaches — oo as C — » C max . 

In this region of large asymmetry it therefore seems more useful to consider the ratio as 
function of the annihilation cross section. This is shown in the right frame of Fig. We have 
chosen a + b/6 as x— axis, so that the cases a = (dashed curves) and b = (solid) can be 
shown on the same scale. The asymmetry C is chosen such that the total relic density has 
the indicated value. We see that the ratio first decreases very quickly when the annihilation 
cross section is increased beyond its minimal possible value, which is reached for C — 0. This 
corresponds to the rapid increase of the required asymmetry C with increasing cross section 
shown in Fig. HJ For cross sections larger than roughly twice the minimal value the curves 
flatten out, and have essentially constant logarithmic slope from then on; this slope depends on 
the desired relic density and on the functional form of the cross section (a— or 6— dominance). 

The product aY x Y x determines the strength of indirect Dark Matter detection signals from 
the annihilation of WIMPs in the haloes of galaxies. In Fig. [6] this product is normalized to 
the same product taken at C — 0, i.e. for symmetric Dark Matter. Notice that the cross 
section in the normalization factor is kept fixed, to the value required to get the indicated 
Dark Matter relic density, while the cross section in the numerator is allowed to vary. Gravity 
acts the same way on particles and anti-particles, and the WIMP distribution in the present 
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Figure 6: Today's annihilation rate for asymmetric Dark Matter divided by the same quantity for symmetric 
Dark Matter (C = 0), as a function of m x ■ C (left) and a + 6/6 (right). Parameters and conventions are as in 
Fig. [5] 



universe is basically only determined by gravityJil The ratio n x /n x should therefore be equal 
to the universally averaged Y x /Y x everywhere. This normalized product thus shows how the 
strength of indirect Dark Matter detection signals varies over the allowed parameter space. 

In the case at hand, even for C = particles can only annihilate with anti-particles. For 
given total Dark Matter density and given annihilation cross section the annihilation rate 
will therefore be two times smaller than for self-conjugate (Majorana) Dark Matter. This is 
essentially compensated by the fact that we need approximately two times larger annihilation 
cross sections for C = to achieve a given total Dark Matter density, compared to the case 
of Majorana dark matter; see the discussion of Fig. |H Fig. can therefore also be used to 
compare with the more frequently studied case of Majorana Dark Matter. 

The left frame of Fig. [6] shows the normalized product a • Y x • Y x as function of m x ■ C, 
which again makes the result almost independent of m x . We also note that, as in the left 
frame of Fig. El the result is almost identical for a = 0, b ^ and for o ^ 0, 6 = 0; this is 
again due to the fact that the annihilation cross section is fixed, for given C, by the desired 
total Dark Matter relic density. 

Recall that an increase of C beyond zero corresponds to an increase of the annihilation 
cross section, albeit initially a very small one (see Fig. Hj). This increases the first factor 
in the product shown here. On the other hand, we just saw in Fig. [5] that this decreases 
the third factor, Y x . Fig. [6] shows that this second effect is always the more important one, 
i.e. introducing an asymmetry in the Dark Matter sector can only reduce the signals from 
Dark Matter annihilation in galactic haloes. However, this reduction is at first very mild. 

^Electromagnetic interactions may affect the distribution of baryonic matter, which in turn can affect the 
Dark Matter distribution through gravitational interactions. 
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For constant decoupling temperature and constant cross section the product Y x ■ Y x scales 
like [y x (C = 0)] 2 — C 2 /4, which has vanishing slope at C — 0. The increase of cross section 
required by the increase of C to keep the total Dark Matter density fixed does not change this 
result significantly: the positive contribution to the slope from the first factor in the product 
is partly canceled by the resulting decrease of the decoupling temperature, which reduces 
Y X {C = 0). Recall also that initially the required cross section depends only very weakly on 
C. 

As a result of this small negative slope at small C, the signals from WIMP annihilation in 
the halo still have about 20% of the strength in the symmetric case when the asymmetry C 
reaches 95% of its maximal value. However, beyond that point the strength of these signals 
plummet quickly, due to the very rapid decrease of the x Te ^ c density which we already saw 
earlier. 

The right frame of Fig. [H] shows the normalized product of annihilation cross section, \ 
density and X density as function of a + 6/6. We now observe an almost constant, negative 
logarithmic slope as the cross section is increased from its minimal value, which corresponds to 
C = 0. Already for small C the essentially vanishing slope in the left frame of Fig. |H] (product 
vs. C) combines with the very large positive slope of Fig. H] (required C vs. cross section) to 
give a finite negative slope. 

In particular, we find that the indirect signals from WIMP annihilation in galactic haloes 
are suppressed by about a factor 10 5 for a cross section about 8.2 (9.0) times the minimal 
allowed cross section, for b = (a = 0). This opens new possibilities for model building. 
A suppression of the annihilation rate by a factor O(10 5 ) below the rate expected for self- 
conjugate thermal dark matter annihilating from an S— wave initial state is required if an 
observed excess of positrons near the galactic center [10] is to be explained [TT] through the 
annihilation of WIMPs with masses near an MeV§ In existing models [UJ this is engineered 
by suppressing annihilation from the S— wave; the velocity dependence of the cross section 
then suppresses today's annihilation rate relative to that in the early universe. Here we see 
that the same suppression can be arranged for asymmetric Dark Matter, if the annihilation 
cross section is about one order of magnitude larger than that required for symmetric DM. 



5 Summary and Conclusions 

In this paper we investigated the abundance of Dark Matter for the asymmetric WIMP sce- 
nario, where Dark Matter particles and anti-particles are distinct. This opens the possibility of 

* Annihilation into e + e~ pairs is only possible if m x > m e . On the other hand, astrophysical observations 
[12] imply m x < 10 MeV in this case. 
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generating an asymmetry between WIMPs and their anti-particles. Here we assume that this 
happened well before the epoch of thermal decoupling of the WIMPs. We do not specify any 
dynamical mechanism for generating this asymmetry; rather, we treat it as a free parameter. 

In our work we have assumed that the Dark Matter has reached thermal equilibrium 
before it decoupled. This allows a fairly generic treatment of the problem. Asymmetric 
Dark Matter is often motivated by explaining the origin of the baryon or lepton asymmetry. 
The baryogenesis or leptogenesis process is usually assumed to happen much earlier than 
WIMP decoupling. We note, however, that there could be cases where the WIMP density did 
not reach full thermal equilibrium. In that case, the abundance evolution would be model- 
dependent. Moreover, we assume that WIMPs can only annihilate with their anti-particles; 
self-annihilation of WIMPs would add additional terms to the Boltzmann equations. 

The Dark Matter decoupled from the rest of the cosmic fluid, and its (co-moving) abun- 
dance froze, when the annihilation rate dropped below the Hubble expansion rate. In the 
asymmetric Dark Matter there are more particles than anti-particles (for the reverse 

case, our results would be the same except that "particle" and "anti-particle" are exchanged), 
the anti-particles annihilated away more efficiently, with large numbers of particles left behind 
without partner to annihilate. As a result, the final abundance is determined not only by the 
cross section as in the symmetric Dark Matter case, but also by the asymmetry. We derived 
approximate analytical expressions for the abundance and decoupling temperature as func- 
tions of cross section and asymmetry, generalizing the well-known treatment of the (chemical) 
decoupling of symmetric DM. 

Numerically we find that the final DM density is already largely determined by the initial 
asymmetry if the annihilation cross section is just two times larger than the value required for 
symmetric thermal DM. We also saw that introducing an asymmetry always reduces today's 
indirect DM detection signals from WIMP annihilation in galactic haloes: even though in- 
troducing an asymmetry requires larger annihilation cross sections to attain the desired total 
DM relic density, the resulting enhancement of the annihilation rate is over-compensated by 
the strong suppression of the final anti-WIMP density. This suppression of today's WIMP 
annihilation rate is desired in models where annihilation of MeV Dark Matter particles into 
e + e~ pairs explains a possible excess of positrons found near the center of our galaxy. 

We find that increasing the annihilation cross section by one order of magnitude suppresses 
today's annihilation rate by about six orders of magnitude. Today's anti-WIMP density 
therefore remains significant only for a rather narrow range of annihilation cross sections 
above the value required for symmetric dark matter. 
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Note Added 

As we were preparing this manuscript, ref. [13J appeared, in which the problem of asymmetric 
WIMP Dark Matter decoupling is also investigated. While the general ideas are similar, the 
main focus of the paper and the details of the treatment are different. In particular, ref. [13] 
presents analytical expressions only for the difference and the ratio of the (co-moving) x an d 
X densities. Numerically our results are in general agreement with theirs. 
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